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GENERAL THEORY OF AIRFOIL SECTIONS HAVING ARBITRARY SHAPE OR 

PRESSURE DISTRIBUTION 

By H. Julian Allen 


SUMMARY 

In this report a theory of thin airfoils of small camber is 
developed which permits either the velocity distribution cor- 
responding to a given airfoil shape, 01 the airfoil shape corre- 
sponding to a given velocity distribution to be calculated. The 
procedures to be employed in these calculations are outlined and 
illustrated with suitable examples. 

INTRODUCTION 

Before the advent of the low-drag and high-critical-speed 
airfoils, the shapes of airfoil sections having desirable aero- 
dynamic characteristics were found by the purely empirical 
method of testing families of related profiles. The pressure 
distribution over any of these shapes could be calculated by 
any of a number of methods, but notably by the method 
of references 1 and 2. 

Subsequently, experimental and theoretical investigations, 
on the one hand, of the laminar boundary layer and the 
phenomenon of transition and, on the other, of the compres- 
sion shock wave promoted a better understanding of the fac- 
tors affecting the drag of airfoils. It became apparent that 
the control of the aerodynamic characteristics of airfoils 
was to be found in the control of the pressure or velocity 
distribution. Hence, in the design of an airfoil having certain 
desirable aerodynamic characteristics, the “inverse” prob- 
lem of finding the shape of an airfoil which would promote a 
specified velocity distribution over its surface became of 
considerable importance. 

One notable method has been advanced (reference 3) for 
solving this inverse problem. However, this method is 
intricate and laborious to employ. 

In this report a new method, which has been used for 
the past several years in the design of a large number of low- 
drag and high -critical -compressibility -speed airfoils, is 
presented. This method, which is comparatively rapid and 
easily applied, may be used to solve either the direct or the 
more important inverse problem. Illustrative examples are 
included. 

THEORY 

It is shown in reference 4 that in a determination of the 
pressure distribution over a cambered airfoil the effects of the 
camber and the thickness distribution may be considered 
independently. Specifically, it is shown that the induced 
velocity at any point on the cambered airfoil may be found 


by superposing the induced velocity at the point due to the 
vortex system, which may be considered to replace the mean 
camber line, and that at the point due to the source-sink 
system, which maybe considered to replace the “base profile.” 
The base profile of the airfoil is the profile if the camber were 
removed and the resulting symmetrical airfoil set at zero 
angle of attack. In the airfoil theory of this report it is 
convenient to consider separately the base profile and the 
mean camber line which together make up a given airfoil. 

CAMBER-UNE THEORY 

Glauert (reference 5) has considered the problem of the 
mean camber line which, in a more convenient form for cal- 
culation and extended so as to include the theory of the strut 
as well as airfoil mean camber lines, is given in the following: 



Figure 1 . — Diagram of mean camber line. 


Consider the mean camber line shown in figure 1. If the 
camber is small, the velocity induced at a point P(x 0 ,y ct) ) on 
the mean camber line by a vortex at any other point P(i,y e ) on 
this line is approximately that which would be induced at 
the point on the a-axis P(z 0 , 0) bv the same vortex at the 

. dr . 

point P(x, 0). If the vortex strength at any point is dx> 

the velocity induced at any point on the camber line due to 
all the vortices distributed along the camber line is 



dx 
' J ’o) 


( 1 ) 


1 



2 


REPORT NO. 833 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


and is perpendicular to the x-axis. The flow direction close 
to the camber line must be parallel to the surface of the 
camber line so that if the angle a between the x-axis and the 
flow direction of the undisturbed stream is small, then 


v dy c 
V 0 ~dx a 


(2) 


where Y 0 is the velocity of the undisturbed stream. 

It is convenient to introduce the new coordinate 0 for x 
such that 


(* \ 
x ~2 (1 — cos 6 ) 1 

x 0 =| ( 1 — COS 0 o) ) ( 3 ) 

c 1 

dx =2 sin 0 dd I 


The coefficients are given by 

A n —-( ^ cos nddd 
ir Jo ax 

The lift force may be found from 

-c P F 0 2 J o ' [4/(l + oos 9)+A 0 "( 1-cos 9) 
A n sin nd sin 0 J dd 
= 7TCpFo 2 (A 0 ' + A 0 "+|a) 


(9) 


where c is the airfoil chord. Assuming the distribution of 
vorticity 0 r (where the prescript 0 indicates that this circu- 
lation applies to an airfoil of zero thickness) along the x-axis is 

2F 0 cot | 9+A 0 " tan ^ 0+ A n sin (4) 

Then 

^ dx=cV o £.Ao' ( 1 -f- cos 9) A- A 0 " (1—cos 0) 


so that the lift coefficient is 

„c,=2»(4,'+^o" +|x.) oo: 

According to theoretical hydrodynamics, in an inviscic 
fluid a strut section with a rounded trailing edge should ex 
perience no net lift at any angle of attack so that in this cas< 

-A„' (11 


+S A n sin n0 sin 0 
i 

Writing sin n9 sin 0=^ [cos (n— 1) 0— cos (n+1) 0 ], then 

from equations (1), (2), (3), and (5) the slope at 0 O may be 
obtained from 



d ?jo_ _1 
dx a T 



(1 + cos 0) +Aq" (1 — c os 0) 
COS 00 — cos 0 


o 2 [cos ( n — 1) 0— cos (n-f- 1) 0] , 

& i 


cos 0 O — COS 0 


dd 


( 6 ) 


whence the relations (9) become 

A ' = l^ i di cosnede 

In the case of an airfoil wherein the trailing edge is sharp 
the “Kutta condition" must be satisfied (i. e., the flow mus 
leave the trailing edge smoothly). To attain this th 
vorticity at the trailing edge must be zero, which require 
that Aq"=0. So the coefficients become 



It is shown in reference 5 that 


C T cos n9 sin nd 0 

I a T“ «0=7T — • — r- 

Ja eos 9 — eos 9 * sin A, 


so that equation (6) becomes for the slope at 0 


(7) 






cos nd d$ 


(13 


dy_c 

dx 


— A„ cos 7i0 

i 


(8) 


It will be noted that for airfoils only A 0 ' and for strut 
only Aq' and Aq" vary with the angle of attack. Th 
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coefficients A n are independent of the angle of attack and 
are functions of the mean camber-line shape only. 

Denote by 6 P the difference at x between the upper and 
lower surface pressure coefficients, Pi—P u . (The pressure 
coefficient is the pressure in terms of g, the stream dynamic 
pressure.) Then, from the Kutta-Joukowski theorem of lift, 


I/ 

fV ° <lx _ 2 rf.r 

q ~V a dx (14 ' 

so that, from equation (4), 

0 P=4 cot \ d+Ao" tan \ 0+S^» sin n ^) (15) 


It has been found convenient in the past to denote that 
part of the chordwise lift distribution which is in magnitude 
independent of the angle of attack and in form dependent 
solely upon the camberline shape, as the basic lift distri- 
bution; and that which is in magnitude variable with the 
angle of attack and in form independent of the mean camber- 
line shape, as the additional lift distribution. (These 
concepts first appeared in reference 6 and were later used 
in the development of the methods of references 7 and 8.) 
Hence, for the infinitesimally thin airfoil or strut the addi- 
tional lift distribution is given by 

0 P a =4 (W cot \ 9+A 0 * tan ~ fl) (16) 

and the basic, by 

0 P ft =42>l* sin nd (17) 


It is convenient to consider the basic lift distribution only 
as characteristic of a given camber-line shape since the 
additional distribution may be modified at will by a change 
in the angle of attack and so, at some angle, must be zero. 
The angle of attack at which the magnitude of the additional 
distribution is zero for an airfoil is called the ideal angle 
(references 6 and 9), and is given by 




1 f'd- 
xjo d 


dye 

dx 


de 


(18) 


The ordinates of the mean camber line corresponding to the 
case when the additional distribution is zero, denoted by y Cbf 
are related to the ordinates y c by 


and so 


From equations 


Vc ^ 

c ~ c c a < 

dy H _dy c 
dx dx a< 

(8) and (17), then 
cos nd 

dx i 



4 


=yVl„ sin nd 

i 


( 10 ) 

(20) 


and the coefficients are given by 

A,=- f ' +? cos nd de=- f ' sin nd de (22) 
n tJo dx IT JO 4 

Using equations (21) and (22), the chordwise list distribu- 
tion corresponding to a given mean camber line* or the mean 
camber line corresponding to a given chordwise lift distribu- 
tion can be found. The calculations will in the general case 
be very lengthy so that it is desirable to replace the Fourier 
expansions by integral expressions, as was done in the de- 
velopment of the method of reference 1. To this end* the 
expression for the Fourier coefficients given by equation 
(22) can be substituted in equations (21). At 0 O then 


and 


Now, 


oPt>o_2 c v<j, ^ g j n cog ^ 
4 x J o dx i 

j-- ~- f ~ X) sin nd cos nd 0 dd 
dx xjo 4 i 

1 


(23) 


sin n0„ cos nd [sin n(0+0o) — sin n(d— 0 o )j 
sin nd cos 7i0 o =^ [sin n(0+0 o ) + sin n(d— 0 O )] 


and further 

n 

y 1 . sin n 


1 /0 + 0\ cos (2»+l) 

^ =5 C<>t ( ~ ) 2si„(*f’-°) 


so that substitution gives 


o Pi 

4 


! 

'=lim 

n->« 

1 

2x . 

nt’KV, 

)-«*(*- 

t) 

l] dd 

c rd y*> 

cos 

(2»+l) (V) 

cos (2 n-\- 

•>(' 


jo dx 


• ( 0+0 o\ 

sin ^ 

0—0 

\ 


sin v 2 ) 

2 

) . 


dd 


and 

d Vc b 


dx 


h I' °x [ cot ( V)+ cot ( VO] 


1 r* o£_* 

2x Jo 4 


cos (2w+l) cos (2n+l) 


sin 


(V°) 


sm 


(VO 


dd 


In the limit, the second integrals in the above relations 
become zero so that the equations may be written as 


A, = j_ 

4 2? r 

d V^_ r 

dx 2 x 


r%Kvo-vvo>) 

j;tKvoMvo>i 


(24) 


and 


( 21 ) 


When either function is known in simple algebraic form, 
it is sometimes convenient to express these integrals as 
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follows: From known trignometric relations, equations (24) 
may be written 

1 C* dVc b sin 9p dd \ 

7T jo dx cos 0— cos 0 O 

(25) 

_i C * oUl sin Odd \ 
t Jo 4 cos 0— cos 0 O ' 


4' 

dx 


which may be useful if the functions under the integrals are 
expressed as simple functions of 0. 

When the functions are expressed in terms of x, the follow- 
ing forms, obtained by substituting the relations of x with 0 
given by equation (3), are sometimes useful: 


dVc, 


<A, 

4 

dye, 

dx 


= _i p <** vx “ (< * o) dx 

xjo (x-Xq )yx(c— x) 


, = 1 f e f P b a 

ir Jo 4(x— 


uPbdx_ 

x 0 ) 


(26) 


The second equation has been used to determine the shape 
of a variety of camber lines, notably the type “a" mean 
camber lines (reference 9) used with the more common low- 
drag and high-critical-compressibility-speed airfoils. 

Unless the algebraic expressions for 0 P b or dyjdx are very 
simple, the direct integrations using equations (25) and (26) 
are not convenient so that, in general, it is desirable to per- 
form the integrations numerically using equations (24) . 

The computation may be shortened considerably by use 
of the following mathematical device: 


£m cot (^±4>) de= - £’f( 2*-e) cot (^°) da 
Hence, equations (24) may be written 

oP» 0 l r 2r dye 


dy c b 

dx 


-ifS- C t 5 )* 

-(tL-ftl. 

(xL 


(27) 


These integrals may be evaluated numerically by the method 
of reference 1 which is given in Appendix A of this report. 

In the preceding theory it was assumed that the airfoil 
was of infinitesimal thickness, hence the velocity at each 
elemental vortex along the camber line was taken to be the 
free-stream velocity V 0 . For airfoils of finite thickness, tbe 
velocity differs somewhat from V 0 . A better approximation 
is to assume that the velocity at each vortex is the velocity 
on the surface of the base profile at the same station. Hence, 


the effect of airfoil thickness will be to change the local lift 
at x to approximately 

p (v0 ' (28) 

where V f is the local velocity on the base profile at x. Tlu 
calculation of V f is considered later in this report. 



Figure 2 .— Additional lift distribution for NACA 0012 airfoil. 



Figure 3.— Additional lift distribution for NACA 0018 airfoil. 


Values of P a lc la calculated by equation (28) for the XAC 
0012 and 0018 airfoils are shown in figures 2 and 3, respe 
tively, along with the values given by the method of reu 
ence 7 which were obtained by interpolation of experimenl 
pressure distributions. Shown dotted are the theoretu 
values pPJcia for the infinitesimally thin airfoil obtain 
from equation (16) and given in table I. In figure 4 t 
calculated and experimental basic lift distributions for t 
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x 
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Figure 4.— Basic lift distribution for NACA 35-215 low-drag airfoil. 

a— 0.5 mean camber line of the NACA 35-215 airfoil are 
shown. It is evident from these figures that equation (28) 
gives a close second approximation to the value of P. 

The basic lift coefficient becomes 

. (29) 

and the quarter-chord moment coefficient is 

w.=J>(£„)(HMD < 3o > 

It is obviously inconsistent to make the approximate 
velocity correction to the lift distribution (equation (28)) 
and not to the velocity ratio v/V 0 in equation (2). How- 
ever, the correction to the lift distribution accounts for 
nearly all the discrepancies between the calculated and ex- 
perimental results so that the additional computational dif- 


ficulties associated with a further improvement of the theory 
are not considered to b.e justified. 

BASE-PROFILE THEORY 4 

The problem of determining the velocity distribution over 
a given base profile or the base profile which will promote a 
given velocity distribution over its surface may be treated 
in a m ann er analogous to that of the mean-camber-line 
theory. 

Consider the base profile shown in figure 5. If the thick- 
ness is small, the velocity induced at a point P (x 0 ,y t 0 ) on 
the surface of the profile by a fluid source or sink at the point 
P(x,0) is approximately that which would be induced at the 
point P(zo,0) by this source or sink. If the source strength 
at a point x is ( dQfdx)dx , then, the velocity induced by all 
sources or sinks distributed along the z-axis will be 



Jo 


dQ 

dx 


dx 


o x 0 — x 


(31) 


The source strength can be related to the shape in the 
following approximate manner: If the profile is thin, the ve- 
locity at the surface does not differ materially from the free- 
stream velocity T 0 , and hence the flow velocity within the 
profile due to the sources and sinks is as a first approxi- 
mation V 0 . Within the profile the difference between the 
quantity of fluid flowing at x+dx and x is the amount sup- 
plied by the source contained within this interval, hence 


so 




dy t 

dx 




dy, 

dx 


and so equation (31) become approximately 



J 


dy t 

dx 


dx 


Xq X 


(32) 


(33) 


Vt 



Replacing x by the 0 coordinate defined by equation (3) 
and assuming that the slope of the profile is given by 

^ = cot ^0+Po" tan \ 0+ jb Bn sin nd (34) 

then by analogy with the similar development in the mean 
camber-line theory 


and 


Yf = B 0 ' — B^' ' — XJ B n cos nd 
Vo i 

fX de 

vjo Vo 

2 r v 

B n = — - 1 -fy cos nd dd 
irj o Vo 


(35) 
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The induced velocity from equation (35) is then 

f=B o '(l + 2cos«)=^0[3-4(f); 


Hence, the ratio of the local to stream velocity is 



ae 

e 


Figure 6.— Theoretical velocity distribution over elliptic base profiles. 



Figure 7. — Theoretical velocity distribution over Joukowski base profiles. 


The velocity distribution for the elliptic base profile i 
calculated by equation (38), and for the Joukowski ba: 
profile as calculated by equation (39) is shown in figures 
and 7, respectively, along with theoretically correct veloeil 
distributions as calculated from potential theory. It is set 
that the approximate velocity distributions are satisfactoi 
except in the region where the slope flyjdx becomes infinit 
This was to be expected siiu'e the assumptions made in tl 
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development of the theory in effect require that dy t jdx be 
small. Inaccuracies due to such infinite slopes may be 
avoided by the following device: It is evident from equations 
(34) and (35) that superposition is permissible in regard to 
both base-profile shapes and the corresponding induced 
velocities. Hence the method may be used to find the 
change in shape corresponding to a specified change in veloc- 
ity distribution or, inversely, the change in velocity distribu- 
tion corresponding to a specified change in shape from some 
“reference base profile.’ ’ If this reference base profile is one 
properly chosen so as to have the same slope characteristics 
at the leading and trailing edges as has the profile under con- 
sideration, then the B 0 ' and B 0 " coefficients in the series to 
represent such changes can be made zero since now slope 
differences need never become infinite. For example, for an 
airfoil with rounded leading edge and pointed trailing edge, 
the Joukowski base profile having the same leading-edge 
radius may be used as the reference base profile; for a strut 
with rounded leading and trailing edges, the elliptic base 
profile having the same leading- and trailing-edge radii may 
be used as the reference base profile. Letting A y t and Av 
represent the change in shape and velocity, respectively, then 
equations (34) and (35) become 


d(Ay t ) 

dx 


m 

sin nd 

i 


B n cos nd 

Vo l 


(40) 


and the coefficients are 


or 



[** d(Ay t ) 
o dx 


sin nddd 


B n = 


2 C* Av , \ 

- rr cos nd dQ 

IT J 0 Vo / 


(41) 


tion. Where the change in shape or velocity distribution is 
known as a relatively simple trigonometric function in d, it 
is sometimes convenient to use the equations , 


and 


d(Ay t< ) = l r* Av sin 9 0 dd 
dx t Jo Vo cos d— cos d 0 

AVo_l C r d(Ay t ) sin ddd 
V 0 ~ir Jo dx cos d — cos d 0 


(43) 


When the change in shape or velocity distribution is 
known as a relatively simple function of x, then it is some- 
times convenient to use the equations * 


A?? 

d(Ay t p) _i p Vp M r x 0 (c-x 0 ) dx 
dx ff Jo (x— x 0 ) ^x(c-x) 


AVq 1 . 

Vo~ T 



d -^dx 

dx 


x—Xo 


(44) 


In the general case when the equations for Av/V 0 and 
d(Ay t )/dx are complex or unknown, the most useful forms 
of the equations are 


defining (^=(£0^ 


Av, 

V t 


rd(Ai/d“] 

"d(Ai/dl 


L dx J r+e 

_ dx J x _j 



f (45) 


These integrals can be evaluated numerically by the 
method given in Appendix A. 


Since Bo' and B 0 " have been set equal to zero, then from 
equation (37) the coefficient B x must also be zero. 

When it is desired to find the change iD the velocity dis- 
tribution corresponding to some given change in shape, the 
given change in shape — if so chosen that d(Ay t )ldx is not 
infinite at the leading and trailing edges — automatically sat- 
isfies the condition that the coefficients B 0 ', B 0 ", and B x are 
zero. However, when it is desired to find the change in base- 
profile shape corresponding to a given change in the velocity 
distribution, the change in velocity distribution must be so 
chosen that 

f T dd = 0 ) 

Jo Vo 

and ; (42) 

C* Am \ 

yr cos ddd— 0 
Jo V 0 I 

if the velocity distribution chosen is to correspond to a real 

base profile. 

As shown in the preceding theory of the mean camber line, 
the sine and cosine series can be replaced by the integral rela- 
tions which are generally superior tor purposes of computa- 


APPLICATIONS OF THE METHOD 

THE CALCULATION OF THE VELOCITY DISTRIBUTION OVER A GIVEN AIR- 
FOIL SECTION 

In this section the procedure to be followed to calculate 
the velocity distribution over a given airfoil section is 
presented, and the calculation of the velocity distribution 
over the NACA 4412 airfoil section is used as an example. 

The procedure may be summarized as follows: The ordi- 
nates of the base profile are obtained from the airfoil ordi- 
nates by removing the camber. The velocity distribution 
over the base profile is found by adding to the known velocity 
distribution over some reference base profile having the same 
leading-edge radius, the change in the velocity dio^ribution 
due to a change in shape from the reference to the given 
base profile. Next, the ordinates of the camber line are 
obtained from the airfoil ordinates by removing the thickness. 
The chordwise lift distribution over this infinitesimally 
thin camber line is calculated and corrected for the effect 
of thickness. Finally, the effect of camber is combined with 
the velocity distribution over the base profile to give the 
velocity distribution over the given airfoil section. 


788741—48 2 
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From the known airfoil section, the ordinates y t c, the 
nose radius of the base profile r L , E Jc, and the ordinates of 
the mean camber line yjc are determined for some or all 
of the standard stations xjc listed in table V . In the general 
case this may be done graphically from a large plot of the 
airfoil, taking care to measure the base-profile ordinates 
perpendicular to the mean camber line at each station. All 
modem NACA conventional and low-drag airfoils are 
formed from specified base profiles and mean camber lines, 
and the ordinates or equations for the ordinates can be found 
in NACA reports. For the NACA 4412 airfoil, the ordinate 
of the base profile (i. e., the NACA 0012) and the equation 
for the ordinates of the camber line are given in reference 10. 

The base-profile velocity distribution is calculated as 
follows: A reference base profile having approximately the 
same nose radius is chosen from the Joukowski profiles listed 
in table II. The thickness ratio of a Joukowski profile 
having the leading-edge radius r L .s. is 

-=0.918 J~ ( 51 ) 

c \ c 

For the NACA 4412, the leading-edge radius is 0.0158c, 
the proper thickness ratio for the Joukowski base profile is 
then 0.1155. It is sufficiently exact and more convenient 
to use the Joukowski section with tjc= 0.12. The difference 
between the ordinates of the given and reference profile is 
found from 


(52) 

c c c 

and listed as in table III. These difference ordinates are 
plotted as a function of xjc, as in figure 8, and the slopes 
graphically determine. These slopes are plotted as func- 
tions of 9 in radians (fig. 9) with the slopes at 9=0 and 9=ir 
arbitrarily set equal to zero. Then the ordinates and slopes 
at the proper 9 stations given in the Appendix A for the 
numerical integration of equation (45) are found and listed 
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Figure 8. — Difference ordinates for base profile of NACA 4412 airfoil section as a functlor 

of z/c. 


in table IV. (Note that the value of the ordinate at statio: 
r+9 must be taken as the value at t— 9 but with opposit 
sign.) If the curve of figure 9 is fair, the “ 20-point” metho' 
of integration, used in the example, is sufficiently exaci 
The integration is performed as illustrated in Appendix A 
The resulting values of Av/Vo are plotted as a function c 
xjc us ing the conversion table V, and the values at the stanc 
ard stations, taken from this curve, are listed as in table II. 
The velocity distribution over the base profile is found froi 


V/__Vr,Av 

vrv 0 +v Q 


(5v 


where the values of V r /V 0 are those for the reference profi 
given in table II. For comparison, these calculated valui 
have been plotted in figure 10 to an expanded scale < 
Vf/Vo along with those determined by the method r< 
reference 1. 



Figure 9. — Slope of the difference ordinates for base profile of NACA 4412 airfoil section as a function of 0. 
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The chordwise lift distribution is calculated, in general, as 
follows: From a curve of the mean-camber-line ordinates as 
a function of xjc, the slopes are determined and plotted as a 
function of 9. From this graph, the ordinates and slopes are 
determined at the proper 9 stations given in Appendix A 
for the numerical integration of equation (27). (The 
values of ordinates at station 0+x must be taken as the 
same in sign and magnitude as the values at station 9— ir.) 
The integration is performed to obtain the value of o P&. 

For the mean camber line of the NACA 4412 airfoil, 
because the equation of the shape is given in elementary- 
functions of xjc, the procedure may be simplified. From the 
leading edge to z/e=0.4 


Vc^L 
c 4 



and from xjc — OA to x/c—l.O 

H[ 0-2 + 0* i-® 1 } 

(It is to be noted that if 0 P b is found from the integral 
equations (25), (26), or (27), regardless of whether y cf> or y e 
is used in the calculation, the additional distribution will not 
appear in the final answer.) Differentiation gives 

"JH D' 8 - 2 (!)} °- 4 <;< lo 


Using equation (3), these may be written in terms of 9 as 
(cos 9-0.2)] O<0<cos~ l (O.2) f 

^=1 (cos 0—0.2); cos _1 (O.2)<0<ir 

* 

These relations could be employed directly in equation (25) 
and the lift distribution obtained. The existence of the 
singular point, however, makes the algebra tedious. 

To employ the numerical method of integration, the slopes 
d/d9 ( dyjdx ), may be obtained by differentiation of the 
above. 


£(&)— S** °<Kcos- ( o.2) 

Te(lfc) = -l sin * cos-(0.2)«K> 

Using the above equations, the ordinates and slopes of the 
curve of dy e /dx as a function of 9 can be calculated directly 
for the proper 9 stations used in the numerical integration 
(table VI) and the integration performed. The values of 
Qp b are then plotted as functions of xjc (fig. 11) and the 
values obtained for the standard stations. These values 
apply to an airfoil of infinitesimal thickness, and the values 
must be corrected to correspond to the airfoil of finite thick- 
ness (loc. cit. equation (28)) by use of the equation (table 
VII) 

P.=(^).A (54) 



Fioure Ll.— Calculated basic lift distribution for NACA 4412 airfoil section. 
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These values are also plotted as in figure 11. The basic or 
ideal list coefficient is then determined from an integration 
of the plot of as a function of x/c. 

The additional distribution can be calculated by finding 
first the values of P a for 0 c <a of unity from 

=!)=(%) (g) (55) 

using the values of (pP a /o c i a ) from table I. This function is 
plotted and the integral c ia ( 0 Ci a = 1) is determined. Then the 
additional distribution is found from 

PA P^l) (M) 

as in table VII. 

The chordwise lift distribution corresponding to a lilt 
coefficient c t can be found from 

P=P»+ (g (57) 



The velocity distribution over the airfoil may be foimd by 
superposition method of reference 4. The upper and lower 
surface velocity distributions, respectively, are 


and 


Yi-Yij p l± 
Vo'-Vo +V,IV 0 

V i Vf P/4 
Fo“F 0 V f /Vj 


(58) 


In table VIII, the velocity-distribution calculations for 
the N AC A 44 1 2 at c { = 0 .72 are given . The calculated values 




of ( V/V 0 ) 2 are plotted in figure 12 along with the calculate( 
values obtained by the method of reference 1 . Also showi 
for comparison in the solid curve are experimental value 
obtained by interpolation of the experimental pressure dis 
tributions of reference 11. 

In figure 13 are shown the calculated values of (V/V Q ) 2 fo 
NACA 4412 at Cj=2.05 as determined by the method of thi 
report and of reference 1. 

The procedure to be followed to calculate the velocity dis 
tribution over a strut section having a rounded trailing edg 
is the same as that for the airfoil, except that the referenc 
profile must.be one having both a rounded leading an< 
trailing edge. The elliptic sections given in table IX ar 
recommended for use as base profiles in these cases. Stru 
sections are usually not cambered, so that the procedure t 
be followed for cambered strut sections is only of academi 
interest and is accordingly not considered here. It maj 
however, be desirable to calculate the velocity distributio 
over a strut section at some angle of attack other than zen 
In this case, the value of 0 P a is calculated from 


2 a (^cot | 


B— tan £ 9 


or 


o P a — 2a 


1 - 2 - 

c 


LVK-DJ 


and the lift distribution from 



(5£ 


(6( 


Figure 12. — Velocity distribution over NACA 4412 airfoil section at ci=0.72. 
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With this lift distribution the integrated value of the lift 
coefficient may not be zero as required by the potential 
theory. A small shift of the distribution would allow this 
requirement to be met, but such a refinement is not justified 
in view of the fact that for real fluids, as a result of the fluid 
viscosity, this requirement is not actually fulfilled. Using 
the value of P a , the velocity distribution over the strut is 
found from equations (58). 

THE CALCULATION OF SHAPE OF AN AIRFOIL SECTION CORRESPONDING 
TO A GIVEN VELOCITY DISTRIBUTION 

In this section the procedure to be followed to calculate 
the airfoil section shape corresponding to a desired velocity 
distribution is presented. To demonstrate the general pro- 
cedure to be employed in the calculation of an airfoil shape 
corresponding to a desired velocity distribution, the shape 
of a “semi-low-drag” airfoil over which the maximum velocity 
occurs near the leading edge on the upper surface and at the 
midchord location on the lower surface is used as a first 
example. As a second example, the specific calculation of 
airfoil shapes having a “double-roof” type velocity distribu- 
tion is considered. 

It should be noted at the outset that the shape of a two- 
dimensional body corresponding to some desired velocity 
distribution may not represent a real airfoil section which is 
both “closed” and pointed at the trailing edge. As a con- 
sequence the desired velocity distribution can be considered 
only as a “first choice” and must be modified, if necessary, 
to satisfy these conditions. The procedure employed in ad- 
justing the desired velocity distribution and calculating the 
shape of the corresponding airfoil is as follows: From the 
desired airfoil velocity distribution, the corresponding ve- 
locity distribution over the base profile is found by averaging 
the upper and lower surface velocities at each chordwise sta- 
tion. This distribution is examined to determine whether 
it corresponds to a real profile and adjusted, if necessary, to 
satisfy this requirement. The base-profile shape corre- 
sponding to this adjusted velocity distribution is calculated. 
The airfoil velocity distribution is finally adjusted, if re- 
quired, to take into account the modifications made to the 
original base-profile velocity distribution. The chordwise 
lift distribution is determined from this adjusted distribu- 
tion. Then the chordwise lift distribution for the airfoil 
with the thickness removed is determined and the mean 
camber-line shape calculated. The calculated mean camber- 
line and base-profile shapes are then combined to give the 
airfoil section shape corresponding to the finally adjusted 
velocity distribution. In the example of the semi-low-drag 
airfoil, these steps are considered in detail. 

Example I 
(a) First choice 

In general, a desired velocity distribution will be one laid 
out to some specified ideal lift coefficient, although it will 
only be required that the quarter-chord moment coefficient 
not be objectionably large. This, as will be seen later, al- 
lows the lift coefficient to be varied within a relatively large 
range without affecting the desired characteristics of the 
velocity distribution. Hence, under these conditions the ad- 


justment to a particular ideal lift may be easily made after 
the base-profile-shape calculation is completed, provided the 
first choice of the velocity distribution is one having a lift 
coefficient within a few tenths of that desired. 

In practically all cases it is required that the shape corre- 
sponding to some desired velocity distribution be one having 
a specified thickness ratio. This requirempnt — together 
with the requirement, discussed previously,, that the desired 
velocity distribution corresponds to that for a real airfoil 
section which is both closed and pointed at the trailing 
edge — complicates the problem since it is not apparent from 
the velocity distribution whether these requirements are 
fulfilled. These complicating difficulties can and '.must 
larger ly be eliminated by choosing the velocity distribution 
wisely. Reference to known velocity distributions over 
existing airfoils having nearly the same thickness ratios and 
similar velocity distributions with that desired will aid in 
this choice. 

Suppose, for example, that the semi-low-drag airfoil, used 
for illustration, is to have a low ideal lift coefficient and a 
maximum thickness equal to 14 percent of the chord, with 
an upper surface velocity distribution similar to that for a 
Joukowski base profile and a lower surface velocity distribu- 
tion similar to that for an NACA 65-series low-drag base 
profile. It is to be expected under these conditions that 
reference to the velocity distribution over a Joukowski base 
profile for which t/c= 0.14, and over an NACA 65-014 base 
profile will aid in the choice of the desired velocity distribu- 
tion. In figure 14, the velocity distributions for these airfoils 



Figure 14.— Preliminary velocity distribution over the semi-low-drag airfoil. 
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are shown. The velocity distribution for the Joukowski sec- 
tion was taken directly from table II of this report. The 
NACA 65,2-016 is listed in reference 9, and the approximate 
velocity for the NACA 65-014, plotted in figure 14, was ob- 
tained using the approximate relation 



Based on the foregoing consideration, a choice is made of 
the desired upper and lower surface velocity distributions, 
designated (' V u jVd)i and (Vi/V 0 )i, respectively, as in figure 
14. " The base-profile velocity distribution (V f /V 0 )i shown in 
the figure is the average of these (loc. cit. equation (58)). The 
subscript i is used to denote that these velocity distributions 
are a “first trial.” 

Having decided upon the desired velocity distribution, the 
shape of the base profile is then determined as follows: A 
reference base profile which has nearly the same leading-edge 
radius as the airfoil to be derived is selected from the Jou- 
kowski base profiles listed in table II. The airfoil to be 
derived in the example will clearly have a leading-edge radius 
approximately midway between that for the Joukowski base 
profile, for which t/c= 0.14, and the NACA 65—014. The 
leading-edge radius of the Joukowski base profile is 

^ = 1.185 =0.02322 

For the NACA 65,2-016 from reference 9 

^•= 0.01704 

c 


so for the NACA 65-014, since the leading-edge radius of any 
airfoil varies as the square of the thickness ratio, 

—^=0.01704 Q4!) 2 =0.013°5 

The airfoil to be found will have approximately the leading- 
edge radius 

r^,_ 0.01305+0.02322 _ 00181 
c 2 

The Joukowski base profile for which t/c= 0.12 has nearly this 
radius (0.01706) and so is used as the reference profile ’in this 
example. 

The difference between the desired and Joukowski base 
profile velocity distributions is found from (table X) 



and the values of (Av/Vo)i cos 0 are calculated using the value 
of the cosines given in the conversion table Y. Then botl 
(Av/V 0 )i and {Av/V 0 )x cos 9 are plotted as functions of 
(fig. 15). 

(b) Adjustment of first choice 
In order that the desired velocity distribution will rep 
resent a real airfoil which closes and has a sharp trailing edgt 
it is required that the relations (loc. cit. equations (42)) 

^d9 = 0 

y o 

f cos Odd — 0 
lo lo 




be satisfied. 
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In the example, numerical integration of the curves of 
fiigure 15 gives 

■ Kfl),— 

J^(pr) cos dde=- 0.0077 


It is clear that if the values of (A»/T 0 )i are decreased over 

the range -^x<0<x, the integrals will both he nearer zero. 

This su g gests that as a second trial, the difference velocity 
distribution (AvjV) 2 shown in figure 15 will more nearly 
satisfy the integral equations (63) and at the same time will 
not destroy the desired characteristics of the given velocity 
distribution. In this case 


J . Ir >=-°- 0034 

cos 0d0= -0.0012 


The conditions of equations (63) are very nearly satisfied. 
The conditions may be satisfied completely by slightly 
translating and rotating the second trial of Av/V 0 . Assuming 
that a small increment 


X$£H' +i2 (H 

be added to the distribution (Av/V Q ) 2 , since 


and 

then making 


and 


f A dd=irki 

Jo \Vo/ 

[a( rr') cos 0 d6=2k 2 

Jo W 




0.0034 

x 


0.0011 


“0.0006 


(64) 


(65) 

i 


the velocity distribution 


(c) Calculation of base-profile ordinates 

The ordinates and slopes of the Av/V 0 curve at the 
proper 0 stations given in Appendix A for the numerical 
integration of equation (45) are found and listed, as* in 
table XI. (Note that the value of the ordinate at station 
tt+6 is equal in sign and magnitude to the value at x— 0.) 
The integration is performed, and the resulting values of 
d(Ay t ) jdx are plotted as functions of xjc (fig. 16) and integrated 



Figure 16.— Slope of the base profile ordinates for the semi-low-drag airfoil. 


to give the values of A y t jc at the desired standard stations. 
The ordinates of the base profile are then found, as in table 
X, from 

Vl=JUjL.*3b (67) 

c c c 


which corresponds to the velocity distribution 


vr 



( 68 ) 


«« 

will completely satisfy the equations (63) as required. 

In table X, the values of 

A 0^=0.0011 + 0.0006 (f-*) 

are given and the difference velocity distribution Av/V 0 is 
calculated, using equation (66). These values are plotted 
as a function of 0 in figure 15. (The value at 0=0 is 
arbitrarily made zero.) 


It is fortuitous that in this example the maximum thick 
ness was precisely 14 percent of the chord as desired. In 
the event that the final thickness were t x and that the desired 
thickness were t 2 , the ordinates and velocity distribution, 
respectively, for the base profile of thickness t 2 could be 
obtained from the equation 



and the approximate equation 




(69) 


( 70 ) 


14 


REPORT NO. 833 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


(d) Final adjustment 

Before the calculation of the shape of the mean camber 
line corresponding to the given lift distribution is undertaken, 
it will be necessary to revise the lift distribution correspond- 
ing to the first trial of the given distribution so as to take 
into account the effect of the changes made to the original 
base-profile velocity distribution to make that velocity 
distribution represent a real profile. This may be done 
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Figure 17.— ' Velocity distribution for semi-low-drag airfoil. 


graphically on the plot of the corrected base-profile velocity 
distribution. 

For example, in the case of the semi-low-drag airfoil used 
for illustration, the upper surface velocity distribution, to be 
similar to the Joukowski type of velocity distribution, must 
be nearly straight from the stations xjc— 0.2 to 1.0, as seen 
in figure 17. The base-profile velocity distribution must be 
a mean between the upper and lower surface velocity dis- 
tributions; hence, when the upper surface distribution has 
been chosen, the lower surface distribution is determined 


(fig. 17). The velocity distribution T over the upper and 
lower surface of the leading-edge section (0<(^/c<C0-2) k 
then suitably chosen so that the base-profile distribution is 
the mean. The basic chordwise lift distribution is related 
to the upper and lower surface velocity distributions by 
(loc. cit. equations (58)). 

Vu Vr, P >/ 4 \ 

v 0 -v^v f /v 0 \ 

F 0 ~F 0 Vf/V 0 1 

(e) Calculation of mean camber-line ordinates 
Now in the calculation of the mean camber-line shape 
(table XII), the basic lift distribution corresponding tc 
zero profile thickness 0 Pb must be used. From equation (54 


o 


P» 


P b 

VflVo 


(72; 


so that this distribution can be obtained directly from 



(73 


To determine the mean camber-line shape which wil 
promote this lift distribution, the procedure is as follows 
The value of 0 P b j 4 is plotted as a function of 0 (fig. 18) 
The ordinates and slopes at the proper 0 stations given ii 
Appendix A for the numerical integration of equation (27 
are found (table XIII) and the integration performed 
(Note that the value of oPb/ 4 at station * -+0 is equal ii 
magnitude but opposite in sign to that at it— 0). Th 
resulting values of dy e Jdx are plotted as a function of xj 
(fig. 19) and the values of y eb obtained by integration. 

The resulting mean camber line is at the ideal angle o 
attack (the angle of attack for which the additional lif 
distribution is zero) and hence, unless the ideal angle i 
zero, the trailing edge is either below or above the xjc axis 
Ordinates of camber lines are generally specified with th 
extremities of the camber line on the xjc axis and designate 
by the conventional symbol y c . 

The ideal angle of attack is simply 


a,=( V -^) -f^) (74 

\ c J xle -o \C J z ; e „i.o 


and 


and 


1/e V'btX 

~ r~ at 


dy c ^dy (j> 
dx dx 


(7i 
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Figure 18.— Quarter basic lift distribution, 0 P»/4 for semi-low-drag airfoil. 



Figure 19.— Slopes of the mean camber line ordinates for the semi-low-drag airfoii. 


In the case of the semi-low-drag airfoil used for the 
example 

<*,= -{- 0.0 1233 radians= 0.706° 

and the values of yjc and dyjdx listed in table XII are 
obtained from the equations 

—=—+0.01233 - 
c c c 


dycjv 1 * 

dx dx 


+0.01233 


Using these values of yjc and dyjdx along with the pre- 
viously determined values of yjc , the ordinates of the airfoii 
can be calculated, as shown in table XIY for the semi-low- 
drag airfoil used as an example, from 


where 


x u x y t . 

— — — — sin 0 
c c c 


&=*-•+*« os/3 
c c c 


*J=*+Vl sini3 
c c c 


!!=£_!/. cos s 
c c c 




(76) 


The resulting airfoil is shown in figure 20. 
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Figure 20. — Semi-low-drag airfoil. 

The chordwise basic lift distribution over the airfoil is 
found from 

Pt= » p *(u) (77) 

and is plotted and integrated to give the basic (or ideal) lift 
coefficient. For this semi-low-drag airfoil, the basic lift dis- 
tribution was plotted and the value of the basic lift coefficient 
was found by integration to be 

c Jft =0.0489 

The velocity distribution over this semi-low-drag airfoil 
was calculated by the method of reference 1 at this value of 
the lift coefficient. The results of this calculation are shown 
in figure 17 to an enlarged scale of V/V 0 for comparison with 
the given (solid line) distribution. It is seen that the 
agreement is close. 

If it is desired to calculate the velocity distribution at a 
value of the lift coefficient other than the basic lift coefficient, 
the procedure to be followed is that given under the heading 
entitled Applications of the Method. 

In the preceding analysis it was tacitly assumed that the 
lift and quarter-chord moment coefficients corresponding to 
the finally adjusted velocity distribution were those desired. 
It is clear that in the final adjustment these coefficients could 
have been adjusted by changing the upper and lower velocity 
distributions, taking care only to keep the average of these 
equal to the base-profile velocity distribution. On the other 
hand, since superposition of camber lines is always permis- 
sible, such adjustments can be made at any time. There are, 
of course, an infinite number of adjustments possible, some 
of which are particularly convenient. As an example, the 
a— l type mean camber line (reference 9) can be conveniently 
used to adjust the lift or quarter-chord moment coefficient 
of this semi-low-drag airfoil. Since the addition of this 
camber line simply shifts the upper and lower surface velocity 
distributions up or down with respect to that of the base pro- 
file, these adjustments do not disturb the desired character- 
istics of the velocity distribution. 

To cite one example, suppose it is desired that the semi- 
low-drag airfoil be adjusted to an ideal lift* coefficient of 0.4. 
With this base profile the a= 1 type mean camber line for 
0 c, 6 of unity attains a lift coefficient c ib = 1.080 (loc. cit. 


equation (29)). The required lift coefficient for this compo- 
nent of the basic lift is 

0.4000 — 0.0489=0.3511 


so that if the mean camber-line ordinates of the a— 1 type 
for oCf b of unity given in reference 9 are multiplied by 


0.3511 

1.080 


0.325 


and the resulting ordinates added to those in table XII, a nei 
mean camber line for which c, 6 =0.4 is obtained. These ii 
turn can be combined with the original ordinates of the bas 
profile to give the corresponding airfoil ordinates. 

As a second example, suppose it is desired that the semi 
low-drag airfoil be adjusted to zero quarter-chord-momen 
coefficient. For the mean camber line given in table XII 
the quarter-chord-moment coefficient is -f-0.0224. The a= 
type mean camber line for 0 c l6 of unity attains a quartei 
chord-moment coefficient of — 0.2506, or, for c me/i ~ 0.022 
the corresponding q c { ft = 0.0894 and Cj 6 =0.0965. Henc* 
zero quarter-chord-moment coefficient can be obtained b 
combining the original mean camber line with an a— 1 typ 
camber line for which , b = 0 .0894. The corresponding basi 
lift coefficient is 


c , =0.0489+0.0965=0.1454 

In retrospect, it can be seen that the more exacting tl 
characteristics of the desired velocity distribution, the moi 
attention must be given to the first choice and final ad jus 
ment of this distribution. In the case of the semi-low-drs 
airfoil used for illustration, it should be quite clear that hg 
the desired lift and moment characteristics both been spec 
fied, the effort required to obtain a satisfactory first, chol 
and final adjusted velocity distribution would be consideiab 
increased. The possible variations in the choice of desin 
velocity distributions are unlimited so that no general rul 
can be laid down for the special treatment required in ea< 
and every case. Facility in the use of this method for t] 
inverse problem can be acquired only through experience. 
Example II 

Experimental studies of a large number of low-drag a 
foils have been made in which the effects of various modific 
tions in pressure distribution were determined. Airfo 
having pressure-distribution characteristics like that of t 
series 3 and 6 low-drag airfoils were found to be definite 
superior in most respects. The somewhat decreasing b 
nearly constant favorable pressure gradients which occ 
over the forward part of such airfoils, from the nose section 
where severe gradients due to the necessary rounded leadi 
edge occur — to the minimum pressure point, are desirab 
This allows relatively large additional lift increments to 
added, as well as some waviness of surface to be tolerat< 
without such additional lift or waviness promoting lo< 
adverse pressure, gradients and so u premature transit) 
to turbulent flow in the favorable gradient region. I 
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nearly constant adverse gradient back of minimum pressure 
has been found to be influential in increasing the critical 
boundary-layer Reynolds number at the transition point 
with' given surface conditions, and so increasing the upper 
limit of the Reynolds number range for lowest drag co- 
efficients. 

Airfoils similar to the NACA series 3, 4, and 6 low-drag 
airfoils are obtained by superposing the ordinates of a base- 
profile shape which promotes a double-roof form of velocity 
distribution on the ordinates of an appropriate Joukowski 
base profile. The problem of finding the shape correspond- 
ing to the double-roof velocity distribution, although it was 
solved originally for the NACA series 3 low-drag airfoils and 
series 4 high-critical-compressibility-speed airfoils by use 
of the numerical method (equation (45) ), is a rather im- 
portant example of one which may be solved by integration 
using the trigonometric expressions for the velocity distri- 
bution in equation (43). In trigonometric form, the equa- 
tions for the double-roof velocity distribution are 


k l +k 2 cos 9 
Vo 

YT—h +&* cos 9 
Vo 


0 <9<d n 

e m <e<r 


(78) 


For those double-roof base profiles the value of Av/V 0 is 
not zero at x=0, so that it is obviously incorrect to super- 
pose such base profiles on- a reference base profile having an 
infinite slope at the leading edge since for such reference 
base profiles V t /Vq=0. Nevertheless, the velocity distri- 
bution calculated by superposition for such combined 
profiles is in reasonably satisfactory agreement, with experi- 
ment except in the immediate vicinity of the leading edge. 

By combining these double-roof base profiles with a 
suitable reference base profile, a variety of satisfactory low- 
drag airfoils can be derived. An example is the super- 
position of a double-roof base profile for which x m =0.4c 
and g 0 4 = 0 . 3059 , and a double-roof base profile in which 
x«=0.7 c and s„. 7 =0.1367 (table XV) on a Joukowski base 
profile for which rf/c=0.10 (table II). A base profile results 
which, when combined with a type a=0.4, c £! ,=0.8201 
mean camber line superposed on a type a=0.7, c {6 = —0.5513 
mean camber line gives an airfoil for which f/c=0.14 and 
which has an upper-surface velocity distribution similar in 
form to that of the NACA 64-series low-drag airfoil and a 
lower-surface velocity distribution similar in form to the 
NACA 67-series low-drag airfoil. This airfoil, the velocity 
distribution for which is shown to an expanded scale of 
V/V 0 by the solid-line curve in figure 21, is completely satis- 


where 0 corresponds the point x and 9 m corresponds to the 
minimum pressure point x m . 

The conditions that the value of Av/V 0 must be the same 
in either equation at 9 m and that the equations (42) be sat- 
isfied require that 

, _£ I" sin 6 m +(r— 9m) cos 0 m ~1 ' 

1 2 L cos 9 m sin 9 m +r—9 m J 

k £ . 

«-2 — 2 

_ . y ( 79 ) 

, s r sin 9 m —9 m cos 9 m I 

3 — 2 [_cos 9 m sin 9 m J rtr—9 m \ 

, _ s cos 9 m sin 9 m —9 m “j 

4 2 |_cos 9 m sin 9 m -\-ir— 0«J < 

where s is the slope of the velocity curve between x/c—0 and 
xjc—x m jc; that is, 




nj 


-(-) 

x/c — x m jc v Ko/ x/c=0 


xjc 


(80) 


The corresponding shape of the base profile which will pro- 
mote this velocity distribution is 


r~ sin 9 m sin 0(1 — cos 9 cos 9 m ) 


c 8 


-(cos 9 m — cos 9) 2 ln 


sin ^ ( 9 m +9 ) 
sin -■ {9 m —9) 


(81) 


cos 9 m sin 0 m +x— 9* 


where the vertical bars indicate the absolute value. In 
table XV, the velocity distribution and ordinates of the 
double-roof base profile which may be superimposed on 
other base profiles are given. 
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Figure 21.— Velocity distribution over a low-drag airfoil having minimum pressure at 0.4 c 
on upper surface and 0.7c on lower surface. 


18 


REPORT NO. 833 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


factory except for the fact that the nose radius is unneces- 
sarily small so that the maximum lift coefficient may be 
adversely affected. One way in which this difficulty may 
be alleviated to some extent is to add the increment A y t 
normal to the surface of the reference base profile rather 
than normal to the z-axis so that the ordinates are given by 

yt=y T + by (82) 

where dy T ldx is the slope of the surface of the reference base 
profile at the station under consideration. This procedure 
is hardly justifiable, however. 

A very satisfactory method for improving the shape at 
the leading edge is to calculate the velocity distribution over 
the base profile carefully by the method of reference 1. 
Then, by using the graphical method of reference 12, a 
change in shape of the base profile and the corresponding 
change in the velocity distribution may be found by trial 
which will allow an increase in the leading-edge radius that 
will not promote a “bump” in the velocity distribution near 
the leading edge. 

DISCUSSION AND CONCLUSIONS 

In the theoretical development of the method of this 
report a number of simplifying assumptions were made in 
order to facilitate the mathematical treatment, some of 
which are clearly contrary to fact. The method of reference 1 
by comparison would appear exact. However, two assump- 
tions common to the development of both methods are 
that the fluid is incompressible and nonviscous. The first 
is justifiable if the velocities are everywhere negligibly small 
in comparison with the velocity of sound. The second can 
never be considered strictly justifiable although in the usual 
Reynolds number range the error is small. It should be 
noted, however, that in the usual case, as may be seen in 
figure 12, the inaccuracy of the method of this report resulting 
from all the other assumptions except that concerning the 
viscosity of the fluid (and compressibility when it is impor- 
tant) is small as compared with the inaccuracies of both 
methods resulting from the neglect of the effects of viscosity 
(and compressibility). 

Concerning the method of reference 1, it has been found 
that the second approximation for the value of e should be 
employed in the calculation if the base profile under con- 
sideration differs markedly from the Joukowski base profile, 


as is the case with a number of low-drag airfoils, particularly 
for thick sections or for those sections with the minimuir 
pressure point far back along the chord. The inadequacy oi 
the first approximation is not very evident in figure 21 
wherein the velocity distribution as calculated by the methoc 
of reference 1 using both the first and second approxima- 
tions is shown for comparison with the method of this report 
since the maximum thickness is fairly far forward and tht 
airfoil is relatively thin. 

The satisfactory application of the methods of this repori 
rests on a thorough understanding of the limitations on th( 
principle of superposition as it applies to the mean cambe: 
line and the base profile. In the theory of the mean cambe] 
line it was assumed that the camber, as well as the slope o 
the camber line, was small. Hence, superposition of mear 
camber lilies or of lift distributions should be permissibfi 
for all usual camber lines, provided the camber or lift i: 
small. Experiment has shown that for usual mean cambe: 
lines calculations based on this method are in good agree 
ment with experiment, provided the basic lift coefficient i 
less than unity, but that even up to basic lift coefficients o 
two the agreement is fair. In the theory of the base profili 
it was assumed, in effect, that the slope of the surface is small 
At the leading edge of an airfbil section and at the leading 
and trailing edges of a strut section, the slope of the profili 
is infinite so that, as was shown for the Joukowski and elliptic 
base profiles, the method of this report cannot be used di 
rectly to determine the velocity distribution or shape of sucl 
bodies. Rather, the method must be used to determine thi 
change in velocity distribution or shape corresponding t< 
some specified change in shape or velocity distribution 
respectively. This change can represent a marked altera 
tion in shape at all points except the leading and trailing 
edges. 

Again, in the theory of the base profile, it was assumed tha 
the profile is thin. Experiments have indicated that th 
method is satisfactory for all airfoils of usual thickness (u] 
to thickness-chord ratios of 0.18) and even reasonably satis 
factory in the case of an airfoil having a thickness-chor< 
ratio of 0.25 (NACA 45-125). 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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NUMERICAL INTEGRATION METHODS 
A numerical evaluation of the integral 

E=l F cot 9 -=* de 

is given in the appendix of reference 1 . 

A “ 20 -point’ ’ solution is 

E=a a (^)+a l (F l -F^)+a 1 (F 2 -F. 1 )+ . . . + 
^(F 9 F _ 9 ) 

wnere F { is the value of F at 


F n is the value of F at 0 o+Jq 

(ti=1, -1, 2,-2 . . . 9,-9) 

(aDo 13 the vaIue of (®) at 9=90 

and the coefficients are 

<*o=0.1000 a 5 =0.0503 

ff!=0.3473 <*8=0.0366 

<* 2 =0.1572 a 7 =0.0281 

<*3=0.0996 <*8=0.0163 

<*4=0.0691 <*„= 0.0080 

The value of &v/V 0 for 0=4 tt/ 1O given in table IV for the 
NACA 4412 base profile, for example, is obtained in the fol- 
lowing cyclic form: 

^r=-|^ 0.1000( + 0.0138) 

+ 0.3473( + 0.0283-0.0192) 
4-0.1572(4-0.0178—0.0006) 

4-0.0996(-0.0044 4-0.0397) 

4- 0.0691 (-0.0406+0) 

+ 0.0503( — 0.0901 — 0.0397) 

+ 0.0366(0+0.0006) 

+ 0.0281( + 0.0901 + 0.0 192) 

+ 0.0163( +0.0406+0.0283) 

+ 0.0080(+0.0044+0.0283)J= —0.0059 

A more accurate “40-point” solution is 

E= 6n ^ ~^~bt(Fi — F-. l ) J rb 2 (F2 — F_ 2 ) + . . . + 6 1( )(.Fi9— F_ l9 ) 


where now 

Fi is the value of F at 


F n is the value of F at 0 O + 

(n=l, -1, 2, -2, . . . 19, -19) 

(jiP) * s va ^ ue 

and the coefficients are given by 

6 „= 0.05000 6 10 = 0.02503 

^=0.34906 6u=0.02139 

6 2 =0. 16129 6 12 =0.01819 

J 3 =0. 10514 6,3=0.01532 

64=0.07735 6,4=0.01273 

65 = 0.06057 6,5=0.01036 

6 fl =0.04918 6, a =0.00814 

6 7 =0,04087 6n=0.00599 

6 8 = 0.03444 6,8=0.00395 

69=0.02929 6,9=0.00197 

The 40-point solution need be employed only when the 
function F changes more or less abruptly with x/c. 
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TABLE I.— CALCULATED VALUES OF ADDITIONAL LIFT 
COEFFICIENT DISTRIBUTION — FOR INFINITESI- 

0 C ‘a 

MALLY THIN AIRFOILS 

o Pa 2(1+ cos d) _2 1 1 — x/c 

0 Ci~ v sin 6 7r"V x 

c * 


X 

~c 

of. 

OCf* 

£ 

C 

0-P . 
*U 

0 

oo 

.4500 

.704 

.0125 

5.658 

.5000 

.637 

.0250 

3.976 

.5500 

.576 

.0500 

2.775 

.6000 

.520 

.0750 

2.235 

.6500 

.467 

.1000 

1.910 

.7000 

.417 

.1500 

1.515 

.7500 

.368 

.2000 

1.273 

.8000 

.318 

.2500 

1.103 

.8500 

.267 

.3000 

.973 

.9000 

.212 

.3500 

. 868 

.9500 

.146 

.4000 

.780 

1.0000 

0 


TABLE II— SURFACE ORDINATES AND VELOCITIES FOR JOUKOWSKI BASE PROFILES 


£ 

C 

</c« 0.02 

(/c— 0.04 

l/c-0.06 

f/e-0.08 

7/C-0.10 

l/c-0.12 

1 

7/C-0.14 

1 

Vr/C 

VrfVt 

Vr/c 

Vr/Vt 

VrlC 

Vr/V* 

Vr/C 

Vr/Vt 

Vr/C 

Vr/Vo 

VrlC 

Vr/V 0 

Sr/C 

Vr/Vo 

0 

.005 

.0075 

.0125 

.025 

.05 

.075 

.10 

.15 

.20 

.25 ‘ 

.30 
.35 
.40 
.45 
.50 
.55 
.60 
.65 
.70 
.75 
.80 
.85 
.90 
.95 
1.00 

0 

.0021603 
.0026358 
.0033772 
.0046853 
.0063737 
.0075003 
.0083121 
.0093443 
.0098525 
,C0 99996 
.0098775 
.0095466 
.0090516 
.0084262 
.0076991 
.0068945 
.0060350 
.0051413 
.0042340 
.0033340 
.0024638 
.0016495 
.0009239 
.0003356 

0 

0 

1.0230 

1.0306 

1.0367 

1.0407 
1.0416 

1.0408 
1.0396 

1.0368 
1.0338 
1.0308 
1.0277 
1.0246 
1.0215 
1.0184 
1.0153 
1.0122 
1.0091 
1.0060 
1.0029 

.9098 

.9968 

.9937 

.9907 

.9876 

.9845 

0 

.0043170 
.0052674 
.0067491 
.0093647 
.0127393 
.0149925 
.0166168 
.0186836 
. 0197030 
.0200000 
.0197586 
.0190994 
.0181107 
.0168610 
.0154070 
.0137980 
.0120784 
.0102902 
.0084744 
.0066730 
.0049314 
.0033015 
.0018491 
.0006716 

0 

0 

1.0045 

1.0321 

1.0556 

1.0732 

1.0799 

1.0799 

1.0783 

1.0736 

1.0676 

1.0616 

1.0564 

1.0492 

1.0429 

1.0366 

1.0304 

1.0241 

1.0179 

1.0117 

1.0055 

QQQA 
. am 

.9933 

.9872 

.9812 

.9752 

.9692 

0 

.0064662 

.0078900 

.0101096 

.0140288 

.0190873 

.0224669 

.0249050 

.0280109 

.0295472 

.0300000 

.0296445 

.0286614 

.0271825 

.0253107 

.0231311 

.0207177 

.0181372 

.0154629 

.0127266 

.0100214 

.0074058 

.0049579 

.0027767 

.0010085 

0 

0 

.9589 
1.0107 
1.0583 
1. 0971 
1. 1145 
1.1171 
1.1159 
1.1096 
1. 1014 
1.0924 
1.0832 
1.0737 
1.0643 
1.0548 
1.0453 
1.0359 
1.0265 
1.0171 
1.0079 
.9987 
.9896 
.9805 
.9715 
.9626 
.9536 

0 

.0086040 

.0104992 

.0134536 

.0186712 

.0254096 

.0299154 

.0331690 

.0373208 

.0393826 

.0400002 

.0395386 

.0382383 

.0362743 

.0337838 

.0308802 

.0276624 

.0242199 

.0206371 

.0169970 

.0133844 

.0098908 

.0066212 

.0037079 

.0013465 

0 

0 

.9016 
.9749 
1.0479 
1.1128 
1. 1452 

1.1522 

1.1523 
1.1454 
1.1350 
1. 1233 
1.1110 
1.0983 
1.0856 
1.0728 
1.0601 
1.0474 
1.0348 
1.0223 
1.0100 

.9977 

.9856 

.9736 

.9617 

.9500 

.9384 

0 

.0107266 

.0130910 

.0167757 

.0232853 

.0316982 

.0373298 

.0414012 

.0466080 

.0492067 

.0500006 

.0494436 

.0478351 

.0453930 

.0422883 

.0386631 

.0346412 

.0303348 

.0258503 

.0212920 

.0167670 

.0123902 

.0082938 

.0046441 

.0016862 

0 

0 

.8427 
.9322 
1.0282 
1. 1210 
1.1719 
1.1849 
1.1872 
1.1807 
1.1685 
1.1542 
1.1388 
1.1230 
1. 1070 
1.0909 
1.0748 
1.0589 
1.0430 
1.0273 
1.0118 
.9965 
.9814 
.9664 
.9517 
.9372 
.9229 

0 

.0128292 

.0156613 

.0200710 

.0278644 

.0379454 

.0447022 

.0495942 

.0558669 

.0590167 

.0600015 

.0593626 

.0574571 

.0545239 

.0508326 

.0464887 

.0416628 

.0364906 

.0311003 

.0256184 

.0201745 

.0149079 

.0099784 

,0055866 

.0020280 

0 

0 

.7875 
.8879 
1.0026 
1.1226 
1. 1946 
l. 2151 
1.2206 
1.2154 
1.2019 
1. 1851 
1.1668 
1. 1478 
1. 1284 
1.1090 
1.0896 
1. 0702 
1.0511 
1.0322 
1.0135 
.9951 
.9769 
.9591 
.9416 
.9243 
.9072 

0 

.0149061 
.0182067 
.0233351 
.0324027 
.0441438 
. 0520252 
.0577410 
.0650927 
.0688103 
.0700034 
.0692988 
.0671100 
. 0637395 
. 0594254 
.0543665 
.0487372 
.0426965 
.0363956 
.0299833 
.0236127 
.0174480 
.0116775 
.0065368 
.0023725 

0 

0 

.7379 
.8449 
.9738 
1. 1190 
1. 2132 I 
1.2428 i 
1. 2524 ! 
1.2495 
1.2351 
1.2161 
1. 1949 
1. 1727 
1.1500 
1.1271 
1. 1043 
1.0816 
1.0590 
1.0368 
1.0150 
.9934 
.9723 
.9516 
.9313 
.9114 
.8919 


—~-0.0(XH7 

c 

^^-0.00190 

C 

^—“0.00427 

c 

TJ ~ -0.00758 
c 

-0.01185 

c 

^^-0.01706 

I C 

■ 

^^-0.02323 

c 


TABLE III.— BASE-PROFILE VELOCITY-DISTRIBUTION 
CALCULATION FOR NACA 4412 AIRFOIL SECTION 


X 

7 

Vj 

c 

V*r 

c 

Ay< 

c 

d(Agt) 

dx 

9 

V, 

Vo 

At> 

Vo 

V, 

Vo 

0 

.0125 

0 

.01894 

0 

.02007 

mm 

— 0. 06t2 

0 

.2241 


-0. 016 

0 

.987 

.025 

.02615 

.02786 


-.0362 

.3176 

1.123 

-.024 


.05 

. 03555 

. 03795 


-.0182 

■llOliM 

1. 195 

-.032 

1. 163 

.075 
. 1 

.04200 

.04470 



.5548 

1.215 

-.034 

1.181 

.04683 

. 04959 


.0012 

.6435 

1.221 

-.033 

1. 188 

. 15 

.05345 

.05587 


.0122 

.7954 

1.215 

-.028 

1.187 

.2 

. 05738 

. 05902 

-.00164 

.0185 

.9273 

1.202 


1. 179 

.25 

. 05941 

.06000 

-. 00059 

.0224 


1. 185 

-.017 

1. 168 

.3 

■BxnEl 

.05936 

.00066 
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TABLE IV.— INTEGRATION CALCULATION OF Av/V 0 FOl 
NACA 4412 BASE PROFILE 
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TABLE IX— SURFACE ORDINATES AND VELOCITIES FOR ELLIPTIC BASE PROFILES 



TABLE X— BASE PROFILE ORDINATE CALCULATION FOR SEMI-LOW-DRAG AIRFOIL 


oooooooooooc 5 o< 


■(CiCC 0 >b*Nt 0 (N 0 » 0 >WS 

5 W©N« 0 ®OMi 0 ^^tDfl >»50 

i^C0-^''3« , oir5»0 , 0-^CCC«*-?QO 

)OOCOOOOOOOOOOO 


>§S8§8©©©©©888 


SggS2882!^SgiS§i^ 


s£S3s3S2§s§sgss§g 

83§6oooSbSSSo8oo8 


gaooooot-wceno-^M-Hoggggg 

sssisisisssslilll 


i i ' f i~\ i i ' ° r r i 


SSS^oS8oS$^S?3o^ 

» i' f u f m r t‘ r 


issifls 


H 8 M 2883 

iwcnSn^®*-* 

IrtH-Upirfciw 


i' f 1 r i' r r * ° i‘ i' f 


38§S© 


iSiSilllll 


§3§§82:S23SaS5 22S 

S 6 b 8 oo 5 oSs$oSooo 


° i i i m t r 


882S£:S2$S8=S$2Sg 

822?5gSS22§gSSS8S 


;|sssgssss|g|||| 


£ 

C 



« 

H 


© 




GENERAL THEORY OF AIRFOIL SECTIONS HAVING ARBITRARY SHAPE OR PRESSURE DISTRIBUTION 
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